Measurement of transport mean free path of light in thin systems 
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We extensively investigate in-plane light diffusion in systems with thickness larger than but comparable with 
the transport mean free path. By exploiting amplified spontaneous emission from dye molecules placed in 
the same holder of the sample, we obtain a directional probe beam precisely aligned to the sample plane. 
By comparing spatial intensity distribution of laterally leaking photons with predictions from random walk 
simulations, we extract accurate values of transport mean free path, opening the way to the investigation of a 
previously inaccessible kind of samples. © 2011 Optical Society of America 
OCIS codes: 290.1990, 310.6860, 290.4210 



The study of light transport in turbid materials is an 
active field of research with many interdisciplinary ap- 
plications in biomedicine [1] , atmospheric study, cosmol- 
ogy [2] , archaeology [3] and also fundamental physics [4] . 
White, diffusive materials are composed of many disor- 
dered particles that perturb the propagation direction 
of light rays with respect to the straight trajectory. The 
usual strategy to measure the transport mean free path 
£, (that is the length after which a light ray loses the 
memory of the incoming direction) consists in compar- 
ing intensity outgoing from the sample with predictions 
of diffusion theory. The most popular techniques are en- 
hanced backscattering cone [5,6] and total transmission 
measurements [7, 8] . The former requires a sample with 
a large scattering volume and the second needs many 
samples with variable volume. Probing in-plane proper- 
ties of thin systems with this two approaches, has to face 
alignment difficulties and reduced signal to noise ratio [9] 
due to losses. £ may be extracted also by measuring spa- 
tial distribution of laterally leaking intensity [10,11], but 
up to now this approach has been limited to large sam- 
ples with i of the order of millimeter. Here we develop a 
measurement protocol that overcomes alignment difficul- 
ties, reduces unwanted scattering produced at interfaces 
and allows to measure £ in systems with thickness of the 
order of tens of micrometers. 

We employ as a test diffusive system a water disper- 
sion of latex particles (Duke Scientific size 520 nm). The 
sample holder (two microscopy coverglasses separated by 
plastic spacers of thickness 5=100 /im) contains a (2 /A) 
drop of diffusive liquid distant about 2 millimeter from a 
drop of 3 /A drop of diethylenc-glycol with 1% volume of 
rhodamine B (see figure lb). As a directional light source 
we used the amplified spontaneous emission (ASE [12]) 
generated between the two coverglasses and in the plane 
of the sample by optically pumping dye doped solution 
with a stripe shaped spot (a 64/im wide and 1.8mm long 
stripe is pumped by a Nd:YAG laser with pulse duration 
of 9 ns, 120 mJ maximum intensity per pulse). The pro- 




Fig. 1. (Color Online) a) image of the light scattered by 
the ASE beam in the rhodamine solution; b) image pf 
two drops inside the sample holder. The shaded rectan- 
gle indicates the pumping area and white line is 870 /im 
long, c) Setup of the experiment. A stripe shaped beam 
from a nanosecond Nd:YAG pulsed laser (PL), shines on 
the rhodamine solution (DYE), producing pencil of light 
emitted at the edges of the inverted area. The ampli- 
fied spontaneous emitted light impinges on the diffusive 
sample (D). Light leaked from the sides of the sample is 
collected on the image plane of a xlO enlarging telescope 
(ET) and collected by the CCD array of a camera. 



duced ASE results in an unpolarised beam characterized 
by a divergence of 2.2° (a picture is reported in figure 
la). Its spectral content is a line narrowed (~ 20 nm 
broad) fluorescence emission centred at 590 nm. 

The beam is aligned to impinge on the drop edge by 
properly orienting the pumping stripe by turning a cylin- 
drical lens (sketch in panel lc). This approach presents 
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Fig. 2. (Color Online) a) Color coded image for the in- 
tensity of laterally escaping light. Graphs b) and c) show 
I(z) relative respectively to the upper and lower shaded 
area of a) . 



several advantage versus direct illumination: I) a perfect 
alignment of the probe beam to the plane of the sample, 
II) a homogeneous illumination of thin sample side III) 
the reduction of artifacts because no optical interfaces 
are present between the diffusive liquid and the light 
source. Moreover contributions from pumping light has 
been eliminated by a colored filter. 

By using an enlarging telescope (xlO magnification, 
0.4 numerical aperture), we project an image of sample 
top, i.e., one of the coverslips, on a CCD camera (Pix- 
elink, pixel size 3.2^m). This image contains the infor- 
mation of how much light is being scattered out of the 
sample at each point (see Fig. 2a). 

Panels 2b and 2c report the intensity / averaged along 
the y direction, over the range Ay~S0fim indicated in 
the figure versus the coordinate z from the two shaded 
area of figure 2a. The two graphs are very similar be- 
cause I(z) is nearly independent of the y coordinates at 
which it is measured if the considered areas lie close to 
the light input position. I(z) presents similar features in 
all samples we measured: a first maximum due to the 
scattering from the physical air liquid interface and a 
second (smoother) maximum due to diffusing photons 
followed by a monotonic decrease. 

In samples much larger than £ the diffuse pattern de- 
pends on the distance between the light source and the 
image plane [11]. In samples in which the thickness is 
comparable with transport mean free path instead, the 
probability for a photon to be ejected at each scattering 
event, that depends on the ratio between £ and S, intro- 
duces losses that affects the shape of I(z). Monte Carlo 
simulations [13], that allows to track photons from the 
source to the detector have been performed to confirm 
this picture. 



In our numerical model random walking particles start 
from a punctual source (located at coordinates [0 0] 
and one mean free path inside the diffusive material [14]) 
and travel performing a isotropic random walk with step 
size equal to £ r until they encounter sample boundary 
(thus escaping from top of the cell) at a certain z. The 
simulation is stopped if either the particle reaches one 
of the lateral edges (located at x=50 and X—-50 /jm) or 
the physical boundary of the input edge (located at z = 
— £ r ) [7] (see figure 3a). The coordinate ze at which the 
random walking particle leaves the sample is recorded. 
A sketch of a single simulation (consisting on a number 
of random walk steps of the order of 10 2 ) reporting the 
trajectory performed by a single photon is represented 
in figure 3a. By repeating the simulation, we obtained 
(in about one day's time of calculation in a standard 
personal computer) 2 • 10 6 values of ze that, organized 
in a normalized histogram, represent an estimate of the 
probability P(z,£ r ) of reaching a certain depth z inside 
de sample. 
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Fig. 3. (Color Online) Panel a) shows the trajectory of a 
single photon particle with a random walk step of R=3 
pjm. On b) P(z) for different values of £ r is reported. 
Graphs are translated vertically an arbitrary amount for 
clarity. 

Some of the curves P{z, £ r ) that have been calculated 
for £ r in the range 1-40 /jm are shown in figure 3b. Each 
simulation included the effect of reflectance at interfaces 
that has been calculated following reference [15] and has 
been taken into account by reducing the escaping prob- 
ability of photons. Although it may be easily taken into 
account [9] , absorption has been neglected because water 
and latex are transparent at the considered wavelength. 
Trying with different kind of sources (gaussian instead of 
punctual), or with more elaborated random walk steps, 
produced differences lower than the simulation noise. 

On the other hand the numerical calculation of 
P{z,£ r ) for a given £ r is performed once and used to 
extrapolate transport mean free path for the considered 



samples simply by rescaling its maximum to coincide 
with the experimental one (Pn(z,£ r ) = aP(z,£ r )). A 
good estimate of £ is the £ r that produces a curve 
Pn(z,£ r ) that best fits the I(z) of the considered sam- 
ple. Figures 4a and 4b show the best fitting Pn(z,£ r ) 
together with I(z) for two different systems: a disper- 
sion of latex beads with concentration of 5% and a pho- 
tonic glass (that is a disordered matrix of latex beads 
blocked in fixed position) obtained by the evaporation 
of the water content from a latex dispersion as described 
in reference [8]. 
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Fig. 4. (Color Online) Measurements and best fits for 
sample "Latex 5%" a) and for sample "PG" b). c) reports 
fit residuals as a function of the spacers size S for the 
sample "Latex 10%". 

Agreement between our experimental protocol and 
standard techniques [8, 16] is confirmed by table 1 in 
which our measurements are reported together with re- 
sults from literature. Samples consist in water dispersion 
of latex beads (520 nm in size) with volume fraction 5% 
and 10%, respectively indicated with "Latex 5%" and 
"Latex 10%". Sample "PG" is the photonic glass. Re- 
ported values and errors result from the statics of the 
measurements on 5 samples for each kind. 

Sample £(/j,m) Measured £(fim) From Refs. 



Latex 5% 


19 ±4 


~21T 


Latex 10% 


8±2 


- 12+ 


PG 


3±1 


-2-3 J 



Table 1. Measured £ compared with literature. T arc from 
reference [16] while " is relative to [8]. 

Graph 4c reports residual of the fitting process (de- 
fined as the absolute value of the area difference between 
the experimental and numerical normalized curves) ob- 
tained measuring sample "Latex 10%" versus the spac- 
ers thickness S. A good agreement is evident from S > 



50/j.m ~6x<, conversely for lower values of S the con- 
tribution from single scattering photons is so strong that 
the shape of I(z) is no longer in agreement with diffusion 
approximations [17]. 

In conclusion by collecting intensity exiting laterally 
in a thin diffusive system shined with light generated 
by in-plane ASE, and by comparing our measurements 
with predictions from numerical random walk simula- 
tion, we are able measure the mean free path. We demon- 
strated how our results are in agreement with measure- 
ments from standard techniques, and that by using our 
approach samples with a thickness up to 6 times the 
mean free path may be characterized. Our measurement 
protocol allows a systematic study of light diffusion in 
intrinsically thin samples like those fabricated by nano- 
lithography or self assembly and optically thin biological 
tissue. 
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